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ABSTRACT 

Using spectroscopy from the Sloan Digital Sky Survey Data Release Seven, we system- 
atically determine the electron density «f and electron temperature of active galaxies 
and star-forming galaxies while mainly focusing on the nan'ow-line regions (NLRs). Herein 
active galaxies refer to composites, low-ionization narrow emission-line regions (LINERs) 
and Seyferts following the Baldwin-Phillips-Terlevich diagram classifications afforded by 
the SDSS data. The plasma diagnostics of and Tg are determined through the /[S ll] 
/I6716//i6731 and /[O ill] /15007//14363 ratios, respectively. By simultaneously determining 
He from [S 11] and from [O III] in our [O III] /14363 emission sample of 15 019 galaxies, we 
find two clear sequences: Tuner k, TcomposHe > Tseyfen > Ts,ar- forming and huner ^ nsfv/«" > 
ncomposite > n star- forming- The typical range of in the NLRs of active galactic nuclei (AGNs) 
is 10" cm . The temperatures in the NLRs range from LO to 2.0x10"* K for Seyferts, and 
the ranges were even higher and wider for LINERs and composites. The transitions of n,, and 
Te from the NLRs to the discs are revealed. 

We also present a comparative study, including stellar mass (M*), specific star formation 
rate (SFR/M*) and plasma diagnostic results. We propose that Yi > Ysy > Yc > Ysf, where 
Y is the characteristic present-day star-formation time-scale. One remarkable feature of the 
Seyferts shown on an M*-SFR/M* diagram, which we call the evolutionary pattern of AGN 
with high ionization potential, is that the strong [O III] /14363 Seyferts distribute uniformly 
with the weak Seyferts, definitely a reverse of the situation for star-forming galaxies. It is a 
natural and well-known consensus that strong [O ill] /14363 emissions in star-forming galax- 
ies imply young stellar populations and thus low stellar masses. However, in the AGN case, 
several strong lines of evidence suggest that some supplementary energy source(s) should be 
responsible for high ionization potential. 

Key words: surveys - galaxies: active - galaxies: ISM - galaxies: Seyfert - galaxies: statistics 



1 INTRODUCTION 

The investigation of the narrow-line regions (NLRs) of active 
galactic nuclei (AGN) is important for a number of reasons. The 
NLR provides a bridge from the familiar and relatively well un- 
derstood to the exotic and poorly known. On one hand, this re- 
gion connects on its inner side with the more compact emitting 
regions that are spatially unresolved (except in the radio) and oc- 
cur in very obscure surroundings. On the other hand, the outer 
boundary of the NLR does not just conjoin with the interstellar 
medium (ISM) of the galaxy b ut it can often be well sp atially re- 
solved at optical wavelengths. iFerland & Netzerl ( 119831) were the 
first to attempt to model the NLRs of AGNs, and many other stud- 
ies followed. Emission-line region models are used to determine 
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the physical conditions o f the ionized gas (for details, refer to 
iDopita & Sutherlandl2003l) . Currently, two main ionization sources 
are most frequently applied: photoionization and shock io nization 
(refer to an overview of modelling NLRs in lGrovesll2007h . In ad- 
dition to the ionization status, the geometry, spatial extent and 
emission from this region could also be explained and reproduced 
by modelling the NLRs. However, before running the models, it 
should be noted that all of these types of models rely on pre- 
estimates of some input parameters (e.g. density and temperature). 
Fortunately, the emission lines observed in the NLRs of AGNs are 
much the same as those observed in H II regions and planetary neb- 
ulae. Therefore, the standard nebular diagnostic methods can be 
applied to determin e the electron density and te mperature in the 
NLRs of AGNs (see lOsterbrock & Ferlandll2006l for more details). 
One of the earliest and most well-studied examples is the NLR in 
Cyg A. The amount of extinction calculated from the Balmer decre- 



©2013 RAS 



2 Z.-T. Zhang, Y.-C. Liang and F. Hammer 



ment giving the best overall fit with the recombination decrement 
forr^=10'' K,w, = 10** cm-3 in Cyg A (.Qsterbrock & Miller 1975.; 
IOsterbrockll98l . lKoskil ( ll978h have found that the average density 
in the line-emitting gas of 20 Seyfert 2 galaxies is approximately 
2000 cm"^ and t hat the a verage temperature is approximately 12 
000 ~ 25 000 K. IWalslj il98j) has studied the Seyfert nuclei of 
NGC 1068 and has found that the electron density range is 9x10* 
to 5x10^ cm"^, assuming an electron temperature of ~15 000 K. 

Two different types of strategies are often applied to investi- 
gate an object class; one is to study one or a few well-observed 
individual objects, and the other is to study a large sample of ob- 
jects. The advantage of the former strategy is the more detailed 
and accurate treatment of data, yet it poses the risk of misidenti- 
fying the truly representative characteristics. In contrast, the lat- 
ter strategy, which is usually based on a homogeneous data-set, 
allows us to utilize the iirformation contained in collective trends 
and correlations. The former strategy has been more frequently 
used when investigating the phy s ical conditions of NLRs (e.g. 



2 SAMPLE SELECTION 



Kraemer. Ruiz & Crenshaw] 1 19981: iKalser et alj |2000| ; ISarth et al.l 



200 ll : ICollins et alj \200^ : (Kraemer et alj |2009|) . In comparison 



there have been only a few previous estimates of the density and 



temperature of a large sample of AGN NLRs (e.g . Rodrigu ez et al 



200C 



Sulentic. Marziani & Dultzin-HacvanI rt et al 

200d IXu et alji2007l) . Because of the success of the Sloan Digi- 



tal Sky Survey (SDSS), we have a large data set to use in order to 
investi ga te the NLRs of AGNs in the local Universe (e.g. Hao et al] 



|2005a" b ; ' Zhou et alj2003 ; iLiu et al.ll2010l : IShen et alj20Ilh . In or- 
der to prepare accurate and trustable input parameters for the next- 
step modelling, we follow the second strategy in this study. 

Narrow emission-line galaxies can be generally separated into 
four classes; star-forming galaxies, composites, low-ionization nu- 
clear emission-line region (LINERs), and Seyfert galaxies. It is 
important to understand how these four classes differ in their 
ionization sources and thus temperatures, densities as well as 
emission, in order to build a pic ture of the local Universe. 
iBaldwin. Phillips & TerlevichI ( Il98lh (hereafter BPT) have pro- 
posed using optical emission-line ratios to classify the dominant 
energy source in emission-line galaxies. The BPT diagrams are sen- 
sitive to the hardness of the ionizing radiation field. The first semi- 
empirical classification to be use d with the standard optical di- 
agn ostic diagrams was derived b y lOsterbrock & de Robertisiri985l 
and lVeilleux & Osterbrocklll987L providing better optical classifi- 
cations. These classification schemes h ave been significantly im- 
proved using the SDSS (e.g. iKewlev et aUi200^ . hereafter Ke06). 
With the improved optical classifications, for the first time, we are 
able to analyse the plasma diagnostic results accurately and to sys- 
tematically distinguish between the different classes within an ob- 
servationally homogeneous sample. 

The present paper is organised as follows. In Section 2, we 
give a description of our sample selection criteria. Sample defini- 
tions are included in Section 3, while Section 4 is devoted to the 
plasma diagnostic results. In Section 5, we focus on NLRs, pro- 
viding statistical and T,, analyses and making comparisons be- 
tween different classes. In Section 6, we further apply the diagnos- 
tic results obtained and present a comparative study dealing with 
the relationships of physical properties, such as stellar mass (M^,) 
and specific star formation rate (SFR/Mj,), in different classes. In 
Section 7, we discuss the shock effects and low-metallicity AGN 
candidates for further study. Finally, we summarize our main re- 
sults in Section 8. Throughout this paper, we assume a flat ACDM 
cosmology with Q„, = 0.3, Q.^ = 0.7, and Ho = 70 km s"' Mpc"' . 



The SDSS jVork et al.ll200ol: IStoughton et" al] |2002|) uses a ded- 
icated, wide-field, 2.5-m telescope ( Abazai ian et al |200 9) at the 
Apache Point observatory. New Mexico, for CCD imaging in five 
broad-bands, ugriz, over 10 000 deg- of high-latitude sky, and for 
the spectroscopy of a million galaxies and 100 000 quasars over this 
same region; these "oals have been realized with the seventh public 
data re lease (DRq Abazajian et al. 2009). The images were re- 
duced jLupton et al.ll200ll ; IStoughton et alj|2002l; fPier et al.ll2003l; 



Ivezic et al 



20041) and c alibrated jHogg et alj 1200 ll ; ISmith et al.1 



2002l ; lTuckeretalJl2006l) . and galaxies were selected in two ways 



for follow-up spectroscopy covering 3800-9200 A at a resolution 
of A/AA =^ 1800. The spectra were taken using 3-arcsec diame- 
ter fibres, positioned as close as possible to the centres of the tar- 
get galaxies. Several improvements have been made to the spec- 
troscopic software since DR5, particularly with regards to wave- 
length calibration, spectrophoto metric calibration and the han- 
dling of strong emis sion lines jAdelman-McCarthv et al] |2008| ; 
lAbazaiian etal]|2009l) . 

The SDSS data base has been explored by several groups, us- 
ing different approaches and techniques. Our analysis is based on 
data products from the catalogue^ obtained by the Max-Planck 
Institute for Astronomy (MPA, Garching) and John Hopkins Uni- 
versity (JHU). The MPA/JHU catalogues are excellent resources 
for statistical studies of nearby narrow-line AGN populations be- 
cause of the very large sample size (e.g. 33 589 AGN in DR2 
and 88 178 in the DR4 catalogue) and because of the inclusion 
of additional data for the host galaxies, such as emission-line 
analyses and some derived phy s ical properties, including stellar 
masses iKa uffman n etal. 2003a; Gallazzi et al. 2005; Sali m et all 
2007|), star formation rates (Brinchmann et al. 2004), and gas phase 
metallicities for star-forming galaxies dTremonti et al.lboojb . The 
MPA/JHU group has publicly released catalogues for a total of 927 
552 SDSS galaxies corresponding to DR7, a significant increase in 
size from their previous DR4 release. A number of improvements 
are made due to develppments both in the SDSS reduction and anal- 
ysis pipelines. The emission-line fluxes given in the MPA/JHU cat- 
alogues have all been corrected for foreground (galactic) reddening 
( lO'Donellll 19941) . 



2.1 Selection criteria 

Our sample selection criteria are as follows. 

(i) The signal-to-noise (S/N) is > 5 in the strong emission lines 
[O II] ^3726, 29 A, H/3 .14861 A, [O III] ^5007 A, [O I] /16300 A, 
Ho- ^6563 A, [N II] ^6584 A, [S II] .16716, 31 A. The S/N crite- 
ria are required so that we can accurately classify the galaxies into 
star-forming or AGN-dominated classes. As discussed on the web- 
sittB the uncertainties (i.e. *_FLUX_ERR) are likely to be underesti- 
mates of the true uncertainties. By using the duplicate observations 
of galaxies to compare the empirical spread in value determinations 
with the random errors, the MPA/JHU group presents the results 
given as 'scale uncertainties' in their DR7 release to correct for the 
underestimates. We adopt this correction and increase the uncer- 
tainty estimates on the emission lines by multiplying by the listed 
line flux uncertainty estimate factors. 



' http://www.sdss.org/dr7 

^ Available at http:// mpa-garching.mpg.de/SDSS/DR7l 

' http://www.mpa-garching.mpg.de/SDSS/DR7/raw_data.html 
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(ii) Objects with no stellar mass (M,.) measurements (both total 
and fibre estimates) and poor specific star formation rate (SFR/M^,) 
estimation (i.e. FLAG = 1) in the catalogues are excluded. We utilize 
these two derived parameters to probe the galaxy formation histo- 
ries of different classes of emission-line galaxies in our s ample. 

(iii) Broad-line AGN contaminations are removed. iHao et all 
l l2005ah have discussed the method of identifying broad-line AGNs 
in the SDSS DR4 spectra. They have proposed FWHM(Hq') > 1200 
km s"' as the selection criterion for defining broad-line AGNs. We 
apply this classification to our data to remove broad-line AGN con- 
taminations, using the SIGMA_BALMER given in the MPA/JHU cat- 
alogues to derive FWHMb„,„,„. = SIGMA_BALMERx 2.355. 

Note that herein we have not included any cuts on redshift, as 
in previous works on SDSS emission-line galaxies (e.g. Ke06). In 
fact, the aperture effects, in Section 5, we present and discuss the 
aperture effects, or more precisely how the determined and 
(plasma diagnostics in Section4) vary as a function of the physi- 
cal aperture size. Although LINERs typically have lower luminosi- 
ties than Seyfert galaxies and are therefore found at lower redshifts 
than Seyferts in the magnitude-limited SDSS, this issue does not 
strongly affect our analysis because we are only interested in the 
typical range and representative values of and T^. 

After applying the above three criteria, our resulting parent 
sample contains 46 867 emission-line galaxies. 



2.2 Extinction correction 

The emission line fluxes are corrected for dust extinction. In this 
study, we use the assumptions given by previous studies. We use 
the Balmer decrement for Case B recombination at T = W* K and 
Hg ~ lO^"'' cm""*, the dust-free Balmer-line ratio Ha/H/? value of 
2.86 for star-forming galaxies and theoretical AGN value of 3. 1 for 
co mposite galaxies, LINERs an d Seyferts, and the relation given 
bv lOsterbrock & Ferlandl (|2006|) 



up lobs 



10 



-cinHa)-f{Hfi)) 



(1) 



Applying the average interstellar ex tinction law , we ha ve f(Ha) - 
f(H/3) = 0.347 as summarized from lFitzpatric3 ( ll999h . 

Theoretically, the Balmer decrement has been determined over 
a range of temperatures and densities for two cases: Case A (as- 
suming the Lyman series is optically thin) and Case B (assuming 
the Lyman series is optically thick). Generally, Case B recombina- 
tion is a better approximation for most nebulae. Some calculated 
Case B results (i.e. ratio descriptions ) for the H I recomb i nation 
lines have be en given by table B.7 i n Dopita & Sutherlan3 ( |2003|) 
or table 4.4 in lOsterbrock & Ferlandl (l2006h . Although we can con- 
versely obtain observational Balmer decrements with Balmer se- 
ries measurements using iig and estimates and the given ratio 
descriptions, the variations as a result of temperature and density 
should not be significant compared with our current values. Fur- 
thermore, although we can predict the HafHfi value for Case B by 
using the and estimates, considerable errors must be intro- 
duced by not only the measurements of line intensity uncertainties 
but also the a lgorithm used for solving and Tg. We should men- 
tion here that lGroves. Brinchmann & Walchej ( l2012l) have found a 
systematic bias in the SDSS DR7 MPA/JPU catalogues (i.e. a 0.35- 
A underestimation of the emission-line equivalent width) and 
have suggested that the DR7 Ay estimates could be overestimated 
by a mean difference of - 0.07 mag. However, this effect is not im- 
portant in our analysis. 



Table 1. Sample divisions based on [O III] ,?4363 
detection quality. Fluxes summarized here are pre- 
extinction correction measurements of the [O III] 
/14363 line (in units of 10"" erg s"' cm"^). Luminosi- 
ties listed are in units of Lq after extinction correction. 
Note that L[oiii] is the luminosity of [O III] A5001. 





Sample A 


Sample B 


Sample C 


Number 


1098 


1409 


13 210 


S/N 


>5 


3-5 


1-3 


Zmedian 


0.070 


0.067 


0.075 


Zrange 


0.03-0.32 


0.02-0.29 


0.02-0.32 


FluX,,ie£^/^„ 


31.0 


13.3 


5.0 


FluX^fl/jge 


6.2-480.9 


4.9-66.1 


1.1-50.1 


log L[oini 


8.32 


8.01 


7.47 


log 


7.59 


7.45 


7.47 



3 SAMPLE DEFINITION 

3.1 Sample: [O III] /i4363 detection quality 

Precise measurements require careful CCD calibrations, along 
with fairly high-quality spectra. In practice, however, the [O III] 
/i5007//i4363 ratio is quite large and is thus difficult to measure 
accurately. Therefore, we further reject 31 150 objects with poor 
[O III] /14363 detection quality (i.e. S/N < 1). The remaining 15 
717 objects are divided into the three samples shown in Table [T] 
strong (sample A), intermediate (sample B), and weak (sample C) 
[O III] /14363 emission samples, respectively. As expected, we see 
that stronger [O III] /I4363 objects display significantly higher lu- 
minosities of [O III] /15007 (L[oiiii). 

The three most widely used BPT diagrams are applied to clas- 
sify the emission-line galaxies: [N II] /i6584/HQ' versus [O III] 
ASmifOp (hereafter [N II]/Hq' versus [O III]/H/3), [S II] /16584/Hq' 
versus [O III] ^5007/1^8 (hereafter [S II]/Hq- versus [O III]/H/3) 
and [O I] ^6584/Htt ver sus [O III] ^5007/H8 (hereaft er [O l\l\la 
versus [O III]/H0) (e.g. [Veilleux & Osterbrock' 'l987|). Using the 
calibrations obtained by |Kewley et al. (2001) (hereafter KeOl), 
iKauffmann et al. I ( l2003cl) (hereafter Ka03) and Ke06, we subdivide 



each of the three samples into four classes: star-forming galaxies, 
composites, LINERs, and Seyferts. 

According to this scheme, our sample of 15 717 galaxies con- 
tains 1842 (1 1.7 per cent) Seyferts, 95 (0.6 per cent) LINERs, 1338 
(8.5 per cent) composites, and 1 1 744 (74.7 per cent) star-forming 
galaxies. The rest are ambiguous galaxies (698; 4.4 per cent), which 
are classified as one type of object in one or two diagrams and 
classified as another type of object in the remaining diagram(s). 
Hereafter we exclude the 698 ambiguous galaxies. Fig.[T]shows the 
sample distributions in the three BPT diagrams. It is confirmed that 
(strong) [O III] /i4363 is more easily detected for Seyferts, which is 
also a clue for some alternative source of ionization in these galax- 
ies. The compositions of the different galaxy classes in sample A, 
B and C are shown in Table. |2] It is clear that Ljoim shows a corre- 
lation with the [O III] /14363 intensity, an effect that is not seen in 
Lh^. Seyferts, LINERs and composites have similar redshift values, 
and are more distant than star-forming galaxies (refer to Table. |2ll. 
In addition, Seyferts and LINERs have higher Lioim than compos- 
ites and star-forming galaxies, which suggests different ionization 
sources are present in different type of galaxies. 

After performing extinction corrections, we calculate Ay in 
seven bins (see Table. [3]l, where Ay = Ry x E{B - V) = 
(3.1c)/1.47 where Ry = 3.1, on average, in the Milky Way and 
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Table 3. Summary of Ay in seven bins. The superscripts sliow tlie sample size of those samples that have a sample size no 
larger than 10. 



logM, (Mo) 


A 


Seyferts 
B 


C 


LINERS 
ABC 


A 


Composites 
B 


C 


Star-forming 
A B 


galaxies 
C 


8.1 






0.29' 










0.25 


0.28 


0.39 


8.7 










0.05^ 




0.53' 


0.29 


0.32 


0.33 


9.3 


0.24' 


0.42^ 


0.35^ 




0.3 1' 


0.24** 


0.37** 


0.32 


0.40 


0.45 


9.8 


0.56 


0.72 


0.47 


1.09^ 


0.57^ 


0.53 


0.60 


0.46 


0.55 


0.68 


10.3 


0.63 


0.62 


0.74 


1.17 


0.57 


0.67 


0.82 


0.68 


0.90 


1.01 


10.7 


0.63 


0.84 


0.84 


1.19^ 0.96 


0.81'" 


0.89 


1.11 




1.08'" 1.36 


11.2 


0.72 


0.86 


0.96 


0.35' 0.95^ 1.10 




1.30^ 


1.35 




0.46 


1.59 
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Figure 1. Three BPT diagrams (from left to right): the [N 1I]/Hq' versus 
[O III]/H8 diagram, the [S II]/Ha versus [O III]/H8 diagram and the [O 
I]/Hq' versus [O IlI]/HyS diagram. Top panels: sample A (blue, upper) and 
sample B (orange, lower) are illustrated. Bottom panels: sample C objects 
(grey) are presented. The KeOl extreme starburst classification hne (solid 
hne), the Ka03 pure star formation line (dashed line) and the Ke06 Seyfert- 
LINER line (dot-dashed line) ai'e used to sepai'ate galaxies into four classes: 
star-fonning galaxies, composites, Seyferts and LINERs. 



Table 2. Number distribution of the subsamples discussed in the text, and a 
summaiy of the mean redshift and luminosity (in units of L©). 



Object Luminosity 
Subsample Number Redshift log L[o iii] log Lh/! 



ASY 


Seyfert 


371 


0.098 


8.64 


7.75 


AL 


LINER 


1 


0.053 


7.66 


7.25 


AC 


composite 


33 


0.107 


8.35 


7.85 


ASF 


star-forming 


529 


0.071 


8.08 


7.44 


BSY 


Seyfert 


419 


0.096 


8.41 


7.55 


BL 


LINER 


5 


0.107 


8.46 


8.01 


BC 


composite 


82 


0.099 


8.07 


7.78 


BSE 


star-forming 


795 


0.073 


7.81 


7.34 


CSY 


Seyfert 


1052 


0.100 


8.24 


7.52 


CL 


LINER 


89 


0.085 


7.74 


7.40 


CC 


composite 


1223 


0.094 


7.79 


7.83 


CSF 


star-foiTning 


10 420 


0.077 


7.37 


7.42 



E(B - V) = c/1.47 fo r the standard Galactic extin ction law (re- 
fer to equation [T] for c: lOsterbrock & Ferlandlliooil) . The follow- 
ing two relations are revealed: comparisons between samples A, 
B and C show Av-a < Ays < Ay-c', comparisons between the 
four classes of emission-line galaxies show Avseyfen < Ay-camp 



< Av-sicirForming S Ay -LINER at coustaut M^,. The first relation could 
be explained by a selection bias because it is far more difficult to 
observe faint lines, such as [O III] /t4363, in a heavily obscured 
medium. Considering the second relation, the sampled LINERs as 
a group have the highest Ay. Because of supernova remnants, the 
high Av in LINERs most likely indicates the coexistence of old 
stellar populations and some specific high-energy process. 

3.2 Class: NLR-dominated objects 

The NLR is generally a f ew kiloparsecs across (e.g. iBennert et al.l 
I2OO2I : iNetzer et al.ll2004h . although ionized gas has been seen at 
even larger ranges in some objects. Therefore, the NLR is a limited- 
small region in its host galaxy, and is most likely a combination of 
gas fluorescing in response to the active nuclei and material be- 
ing ionized by massive stars nearby. The SDSS 3-arcsec diameter 
fibres subtend a variable amount of the total light as the redshift 
varies; therefore, the aperture effects are significant when consider- 
ing physical conditions such as n,, and T^. In this section, we focus 
on the NLRs in our samples to gain insight into the relative contri- 
butions of the NLR and disc as the physical aperture sizes vary. 

In the spectra of Seyfert 2 galaxies, the narrow (permitted 
and forbidden) emiss ion lines sho w typical linewidths (FWHM) of 
300-500 km s"' (e.g. lKonatschiiv & Wans 200^). For LINERs, the 
FWHM (permitted and forbidden) could be even greater than those 
of the Seyfert 2 galaxies. Fig.|2]shows histograms of FWHMB„;,„e, 
for active galaxies (i.e. Seyferts, LINERs and composites; cf. the 
estimation of FWHMBn/,,,^,- in Section 2). The linewidths range from 
100 to 600-800 km s"' with median values of 250, 335, and 419 
km s"' for composites, Seyferts, and LINERs, respectively. We 
also check the FWHM of forbidden lines and the median values 
(in the same order) are: 264, 341 and 439 km s"'. Not surpris- 
ingly, the transition objects (i.e. composites) have narrower lines 
compared to LINERs, as is expected because of differences in their 
average Hubble types and the well- known de pendence of nebular 
linewidth on bulge prominence (e.g. lHciiri996l) . In addition, we see 
that LINERs have wider permitted lines than Seyferts. Table|4]sum- 
marizes the FWHM-selected objects in detail. It is not surprising 
that objects with V'^HMsaimey > 300 km s"' are much more easily 
found in LINERs and Seyferts than in composites. Here, we no- 
tice that there are 168 (with S/N > 1; ~1.4 per cent) narrow-line 
star-forming galaxies. Thus we plot these star-forming galaxies in 
the BPT diagram (i.e. [N Ill/Ho- versus [O III]/IIy8) and we find that 
these objects generally lie along the Ka03 pure star-forming line 
with no significant dispersions; we suppose they could be misclas- 
sifie d and actually b e tran sition objects (i.e. composites). 

iBennert et al] 1 I2OO6I) have studied the NLRs in Seyfert 2 
galaxies from spatially resolved optical spectroscopy. They have 
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Figure 2. Histograms of FWHM of Balmer lines for the thi'ee types of 
emission-line galaxies (from top to bottom): composites, Seyferts and LIN- 
ERs. 



Table 4. Summary of objects with FWHMb„/,„„ > 300 km s"'. 




Object 



S/N > 1 S/N > 3 

Number Per cent Number Per cent 



Figure 3. This ternary diagram illustrates selection criteria of the NLR- 
dominated objects (top red diamond). The disc-contaminated NLR objects 
occupy the left (L; green triangle) and right (R; orange trapezoidal) regions, 
while bottom blue trapezoid is occupied by the non-NLR objects. The ar- 
rows show scaling directions. Refer to equation (|2) for the power law index 



Seyfert(1842) 


1141 


61.9 


540 


29.3 


LINER(95) 


82 


86.3 


6 


0.63 


Composite(1338) 


383 


28.6 


46 


3.4 


SF(11 744) 


168 


1.4 


28 


0.2 



suggested that the lower limits of NLR size are 1-5 kpc and are lim- 
ited by either the S/N data or the lack of a strong surrounding stel- 
lar ionization field. We convert redshifts int o physical aperture size 
(i.e. (j) in kpc) for the SDSS 3-arcsec fibre. IHoI ( l2008h has claimed 
that the NLRs in lu minous Seyfer t s span from 50 pc to 1 kpc in 
radius. Additionally. iBennert et al. I( l2002h have measured the NLR 
sizes based on a sample of Seyferts. They have found a strong cor- 
relation between the NLR radius (Rnlr) and t he [O III]/t5007 lumi - 
nosity, which was scaled to the cosmology bv lNetzeret"al] ( l2004h . 
adopted in this paper. This can be written as 

Rr^LR = 2AL°i^,^^(mkpc), (2) 

where L^q^j^IJ,* = Lfoim / lO"*^ erg s"' (the uncertainty on the con- 
stant 2.1 kpc is of the order of 15 per cent). We adopt this correla- 
tion to constrain the upper limits of Rnlr in order to compare with 
the sampled aperture size (Robs)- Fig. [3] illustrates our NLR selec- 
tion criteria. We require that ln0 < 1 kpc, FWHMg„;,„j,- > 300 km 
s"' and Robs < Rnlr for the so-called NLR-dominated (hereafter 
ND) objects. The other two classes are the disc-contaminated NLR 
(hereafter DN) and non-NLR (hereafter NN) objects. Table |5]is a 
summary of the three classes of objects. We examine how and 
Tg vary as velocity dispersion increases. T he results agree with the 
suggestion given bv lHeckmanetalJjl98lh that the kinematics and 
physical conditions of the narrow-line emitting gas are not strongly 
coupled. We utilize the class definitions in Section 5. 

4 PLASMA DIAGNOSTIC RESULTS 

The two best examples frequently used to measure the electron den- 
sity are [S II] A6716, 31 and [O II] A3126, 29. However, the typical 



linewidths (FWHM) in NLR s, ranging from 300 to 500 km s"' (e.g. 
iKollatschnv & Wansl[200^) . are comparable to, or are even larger 
than, the separation of the two [O II] lines, approximately 300 km 
s"'. Therefore, the [O II] intensity ratio, which is a good electron 
density diagnosti c in H II regions and planetar y nebulae, cannot be 
applied to AGNs l lOsterbrock & Ferlandll2006l) . The [S II] ratio has 
two practical advantages for the determination of n^: for z < 0.3 the 
two [S II] lines fall in the rest-frame optical range and are hence eas- 
ily observed, and the lines have a small enough separation and are 
thus not sensitive to reddening correction. The [O III] /^5007//^4363 
ratio is sensitive to density as well as temperature. The 5007 and 
4363 lines have different critical densities, and at densities of > 10* 
cm"\ the 5007 line is heavily supp ressed because of the colli - 
sional de-excitation of the 'D2 level dOsterbrock & Milleij[l975l) . 
Because of the different critical densities of the line transitions, the 
[O III] ratio is no longer a valid thermometer in the high-density 
(i.e . >10* cm ~^) regions usually seen in broad-line regions (BLRs; 
e.g. lHcil2008l) . However, considering that the average physical aper- 
ture coverage of the SDSS 3-arcsec diameter fibres for galaxies at z 
> 0.02 (cf. Table[T) is much larger than the typical size scale of the 
BLRs, the [O III] ratio is therefore still a reasonable temperature 
estimator in our case. Moreover, although the [O III] ratio is rela- 
tively sensitive to reddening correction because of the significant 
line separation, the mean impact for temperature measurements is 
no more than 10 per cent in our sample. In this study we choose 
to derive the electron density and electron temperature using I[S II] 
.i6716M6731 (hereafter R[S II]) and /[O III] .15007M4363 (here- 
after R[0 III]), respectively. 

The calculations of electron density and temperature are pro- 
cessed with the current atomic dat^l using the TEMDEN task 



References are listed at |http://www.noao.edu/noao/staff/shaw/nebul"M] 
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Table 5. Mean values of redshift, Ay, FWHMBaimer (FWHMg; km s ') and luminosity (Lq) of the three classes: NLR- 
dominated (ND), disc-contaminated NLR (DNL and DNR; refer to Fig.|3) and non-NLR (NN) objects. 



Class 


Number 


Redshift 


Av 


FWHMb 


log L[o nil 


log ill/, 


Sample A 


Sample B 


Sample C 




/J 


U.Uj J 


1 .UJ 












9n 
zu 


NDt 


2 


0.033 


2.16 


497 


8.69 


8.09 






2 


NDc 


21 


0.034 


1.67 


379 


8.09 


8.17 




2 


19 


NDsF 


1 


0.041 


2.20 


306 


8.21 


8.86 






1 






118 


w.oo 


S94 


8 34 


7 54 


104 


1 




DNLl 


429 


0.096 


0.94 


451 


7.82 


7.50 


1 


3 


65 


DNLc 


69 


0.118 


0.97 


378 


8.06 


8.04 


9 


30 


274 


DNLsF 


11 


0.140 


1.20 


351 


8.30 


8.47 


3 


13 


118 


DNRsY 


313 


0.108 


1.11 


449 


9.00 


8.16 


139 


107 


183 


DNRz. 


49 


0.076 


1.66 


494 


8.66 


8.46 




2 


9 


DNRc 


134 


0.102 


1.75 


373 


9.00 


8.95 


2 


3 


44 


DNR5/7 


33 


0.122 


1.39 


335 


9.01 


8.97 


7 


51 


21 


NNsy 


701 


0.084 


0.66 


234 


7.95 


7.24 


98 


152 


451 


NNl 


13 


0.048 


0.84 


234 


7.01 


6.71 






13 


NNc 


955 


0.087 


0.90 


218 


7.68 


7.70 


22 


47 


886 


NNsf 


11 576 


0.075 


0.68 


151 


7.45 


7.41 


519 


111 


10 280 



within the IRAP|1 STSDAS packagfH jShaw & Dufou3l 19951) . The 
task is based on the five-level ato m program first developed by 
Ide Robertis. Dufour & Huri3 ( Il987h . including diagnostics from a 
greater set of ions and emission lines, particularly those in the ultra- 
violet satellite (e.g. International Ultraviolet Explorer and Hubble 
Space Telescope) archives. 

We treat samples A, B and C in a uniform way, except for sam- 
ple C, which had poor [O III] /14363 detections. That is, first we 
apply the uncertainty measurements of [O III] /14363 to calculate 
the upper limits of the R[0 111] ratio at a Icr uncertainty, before per- 
forming the plasma diagnostics. In other words, the sample C tem- 
peratures are correspondingly lower limits because higher R[0 111] 
ratios lead to lower T,,. Higher assumed temperatures correspond 
to lower electron densities, and vice versa. Our plasma diagnostics 
are calculated as follows. 

(i) The initial estimate of n^ is obtained with the [S 11] ratio as- 
suming, for example, r=10 000 K (for those with 1.43 < R[S 11] < 
1.47, some lower values of T are assumed). 

(ii) This density is then used to derive T^,[0 111]. If this differs 
from the initial value of T by more than 10 per cent, a new value of 
«(, is then derived using the T^[0 111] obtained. 

Note that because the density-sensitive tracer R[S II] used here is 
easily excited in a range of physical environments, we suggest that 
the ;ie estimates should be considered as the mean values of the 
sampled regions because of the multi-phase nature of the ISM. Cor- 
respondingly, the value of derived from ^[O 111] with n^ from 
i?[S II] will be biased to higher values. Given n,, < 10^ cm"', the 
ofi'sets wou ld generally be no more than 2000 K as long as < 20 
000 K (e.g. lOsterbrock & Ferlandll2006l) . Therefore, we claim that 
the diagnostics applied compose the most acceptable and feasible 
methodology for analysing a large SDSS sample such as ours. 

At an assumed temperature or density, there are some objects 

' IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in As- 
tronomy under cooperative agreement with the National Science Founda- 
tion. 

* Available at http://www.stsci.edu/institute/software_hardware/stsdas 



Table 6. Summary of objects that have ' T,, with [O III] and Ng with [S 
II]' (TONS). The four cases are: TONS], the calculation using the two 
hne ratios leads to no consistent result; TONS2, R[S 11]>1.47, and thus 
there are no density results; TONS3, R[0 11I]<11.1, and thus there are 
no temperature results; TONS4, both density and temperature determi- 
nations fail. 





TONS] 


TONS2 


TONS3 


TONS4 


Total 


ASY(371) 


12 


1 


1 




14 


AL(1) 












AC(33) 


3 




2 




5 


ASF(529) 


96 


47 






143 


Total(934) 


111 


48 


3 




162 


BSY(419) 


17 


5 


1 




23 


BL(5) 




1 






1 


BC(82) 


14 


3 


5 




22 


BSF(795) 


204 


104 


14 


2 


324 


Total(1301) 


235 


113 


20 


370 




CSY(1052) 


51 


37 






88 


CL(89) 


13 


2 


2 




17 


CC(1223) 


181 


71 


17 


1 


270 


CSF(10 420) 


3462 


2328 


101 


35 


5926 


Total(12 784) 


3657 


2438 


120 


36 


6251 



labelled as 'TONS', which is short for with [O 111] and N, 
with [S 11]' ; this designation means that no (consistent) results were 
found using the simultaneous calculation processes (see Table |6]l. 
For example, it is noticeable that 2541 CSF (i.e. sample C star- 
forming galaxies) show > 1.5x10'* K (cf. columns 1 1 and 12 of 
Table |7j, while 1742 (~ 70 per cent) of these are actually TONSi 
(1109) or TONS2 (633) objects. Therefore, we further check the 
intensity ratios of these TONS objects. 

For the TONS2, TONS3 and TONS4 objects, using the uncer- 
tainty measurements of the [S 11] and [O III] lines to calculate the 
lower limits of R\S II] and the upper limits of R[0 111], we find that 
~ 80 per cent of both TONS2 and TONS3 along with all TONS4 
objects show computable ratios. For the remaining 20 per cent of 
TONS2 and TONS3 objects, we suggest that the uncertainties in- 
troduced by [O III] extinction corrections have the greatest impact. 
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Table 7. Summary of the physical conditions (density and temperature). Note that by applying our determination procedures, the sample C temperatures are 
the lower limits. Ns denotes the number of objects with R[S II] > 1.47 (i.e. the upper limit for TEMDEN calculation using i?[S II]). No denotes the number of 
objects withR[0 III] < 11.1 (i.e. the lower limit of R[0 III] for TEMDEN calculation). j (per cent) is the number of objects with 1.5 < < 2.0 (10* K) and 
the fraction of the coiTesponding subsamples. A'2.o(per cent) is the number of objects with Tg > 2.0 (10'* K) and the fraction of the corresponding subsamples. 
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Figure 4. NonnaUzed histograms of electron density (a) and electron temperature (b). For each panel, from top to bottom, we show Seyfeils (red solid) and 
LINERs (orange dot-dashed), composites (darkgreen), and star-forming galaxies (grey), and from left to right, we show samples A, B and C. 



Therefore, these three cases are most liliely introduced by observa- 
tional uncertainties. In the TONSi cases, we suggest an additional 
possibility, that the intrinsically different spatial distributions of the 
[S n] and [O III] emitting regions might lead to inconsistent results 
in the calculations of and T^. 



The plasma diagnostic results are summarized in Table |7] and 
Fig. |4] shows the normalized histograms of and T^. The main 
results are as follows. 

Electron density. Considering the density results, we can see a 
clear sequence (cf. columns 5 and 6 of Table |7}: nuNER <: 'i-Seyferi 
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> ficomposiie > nsiar- forming ■ The typical 11^ Uncertainties are 150 cm"' 
for Seyferts, LINERs and composites, and 80 cm"' for star-forming 
galaxies. Comparing samples A, B and C, the CSY and CL pos- 
sess significantly lower (~150 cm"') than A/BSY and A/BL, 
most likely indicating some effect of shocks in strong [O III] /i4363 
emitters (please refer to Table|2]for the abbreviations). The charac- 
teristic range of active galaxies, including Seyferts, LINERs and 
composites, is 10^"' cm"', while that of the star-forming galaxies 
is 10'"- cm"'. As shown in Fig.|4ja), the sequence mentioned 
above is visible. However, we see a bimodal number distribution 
of which is especially significant for star-forming galaxies. To 
check whether this is a real effect, we examine the ^[S II] distri- 
bution and we can confirm no corresponding bimodality. In fact, 
the typical uncertainty for objects with rig < 10 is ~40 cm"' in our 
sample. As a result, the bimodal distribution is most likely intro- 
duced by the algorithm used for calculating the electron density in 
the TEMDEN package and is therefore not a physical real effect. 
Electron temperature. The temperature results are listed in Table 
|7](cf. columns 9 to 12). Fig. |4jb) shows normalized temperature 
histograms for each galaxy type. As given above, sample C tem- 
peratures should be viewed as the lower limits for comparison. 



First, a sequence is clearly shown: Tl_ 



> T 



> T: 



Seyfer! 



> Tsiar-fonniiig- The typlcal uncertainties of for the four galaxy 
classes are (in units of lO** K) 0.2 for Seyferts, 0.3 for LINERs, 
0.5 for composites and 0. 1 for star-forming galaxies. Secondly, the 
sampled LINERs display a mean of approximately 2x10'* K (cf. 
column 10 of Table|7) and show the highest number fraction of high 
Tg species (i.e. > 1.5x10" K; cf. column II and 12 of Table[7ll. 
Thirdly, a remarkable result is that the composites as a group show 
a value of Tg very close to that of the LINERs and approximately 
5000 K higher than that of the Seyferts. This result reveals the ex- 
istence of very different dominating ionizing source(s) in compos- 
ites from those in star-forming galaxies. For star-forming galaxies, 
the derived temperatures are generally consistent with the classical 
diagnostic results of bright planetary nebula (~ l.lxlO'* K). More- 
over, as revealed in Fig. |4] the active galaxies (e.g. LINERs and 
composites) show a significant tail extending to much higher tem- 
peratures, up to 6x 10"* K, which is far too high to be explained only 
by photoionization caused by stars. 



5 ELECTRON DENSITY AND ELECTRON 
TEMPERATURE IN NLRS 

Electron density in NLRs. As shown in Fig. [5] most NLRs lie in 
the range of 10^"' cm"'. Only seven Seyferts, seven composites 
and one LINER have lower values of n^, some of which might be 
the result of flux measurement uncertainties. We have tried higher 
densities of 10^"'' cm"' to calculate Tg[0 III] for our selected NLRs, 
but a significant fraction of Seyferts show values no greater than 
6000 K, which could not be repres entative of the cases of NLRs of 
AGNs found in the literature (e.g. iBennert et al.ll2006l) . Thus, our 
results suggest that the characteristic range in NLRs of active 
galaxies is 10^"' cm"'. 

Electron temperature in NLRs. As illustrated in Fig. |5] the typi- 
cal range of Tf, in Seyferts is 1.0-2.0x10'* K, with a median value of 
1.4x10" K. Because both the selected LINERs and composites con- 
tain larger fractions of sample C objects than the Seyferts, and the 
Tg derived could thus be biased to lower values, the typical range 
of these two classes could be higher and wider than that shown in 
Fig.H] 

Fig-IHplots ng (top) and Tg (bottom) as functions of physical 



o 



~1 1 [ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 T" 



Seyfert 

LINER 
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I 6f 




log (cm ^) 

Figure 5. Electron density versus electron temperature. Sample distribu- 
tion of the NLR-dominated objects: Seyferts (blue boxes), LINERs (red 
trianglesjand composites (green circles), with error bars showing Icr dis- 
persions. Open symbols denote objects from sample C with lower Tg limits. 



aperture size ln0 for Seyferts in bins of ln0. For a precise evalua- 
tion, we further calculate the mean values of ng and Tg for different 
types of objects. As shown in Table [8] the ng and Tg of Seyferts 
show transitions from the NLR to the disc: n^osy > noNRsy > ^dnlsy 
> nNNsy and TNOsy > ToNRsy > ToNLsy > T^Nsy (refer to Table |5] 
for the classifications). The unexpected low value of Tg for NLR- 
dominated NLR LINERs and the high value of Tg for non-NLR 
LINERs can most likely be explained by selection biases: the for- 
mer is a result of the small sample size (i.e. only two from sample 
C) and the latter is because the non-NLR LINERs are sampled at 
lower redshifts than the disc-contaminated NLR LINERs and are 
biased to regions close to the galaxy centres. For composites and 
star-forming galaxies, we find no clear transition, as seen in the 
Seyferts. Therefore, we have checked further and found that this 
result was due to misclassificat i on; tha t is, the correlation (see equa- 
tion|2j found bv lBennert et al.l ( |2002|) was derived from Seyferts. 

We have reported single and statistical estimates of ng 
and Tg for the NLRs of active galaxies. Physically, the 
emission-line regions are, by their nature, multi-phase me- 
dia, and thus densities and temperatures are strongly coupled 
with non-uniformity. When modelling the NLRs, both den- 
sity and temperature could vary with different physical con- 
ditions, and such n on-uniformity or multi-phase nature should 



be considered (e.g. Binette. Wil son & Sto rchi-Bergmanril 1 1 9961 : 



iKomossa & Schulzlll997l : lBaskin & Laor 200^. We will leave this 
discussion to a companion paper, which will focus on modelling. 



5.1 High-r^ NLRs in line-ratio diagnostic diagrams 

One item remaining in the problem of modelling of NLR is the so- 
called temperature problem: the difficulty in re producing high Tg[0 
III] w ith photoionization models. For instance, jHeckman & Balickl 
[197^ have claimed that Tg > 20 000 K requires a source of en- 
ergy in addition to photoionization. Traditionally, two different av- 
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Figure 6. Physical aperture size ln0 versus itg (top) and (bottom) for 
Seyferts in bins of ln0. Empty boxes are the median values oflrup, and T^, 
with error bars showing lo" dispersions: NLR-dominated (ND, red), disc- 
contaminated NLR (DNR, orange; DNL, green thin; refer to Fig. |3) and 
non-NLR objects (NN, dashed blue). 



Table 8. Mean jv and of the NLR-dominated (ND), the disc- 
contaminated NLR (DNL and DNR) and the non-NLR (NN) objects. 
Strong (sample A), intermediate (sample B) and weak (sample C) [O 
Ill]/i4363 objects are all included. For classifications, please refer to Fig. 
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ND 


DNR 


DNL 


NN 


ND 


DNR 


DNL 


NN 


SY 


427 


394 


306 


255 


1.40 


1.39 


1.28 


1.29 


L 
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244 


92 


1.07 


1.75 


1.73 


1.87 


C 


197 


166 


151 
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1.54 


1.63 


1.71 


1.61 


SF 


109 


148 


94 


44 


1.78 


1.44 


1.93 


1.32 



enues have been explored when modelling the NLRs: photoioniza- 
tion and shock ionization. It is generally believed that photoioniza- 
tion is the dominant excitation mechanism i n most AGNs, wh ile 
the intrinsic defect of all shock models (e.g. lAllen et al]|2008h is 
that they require shocks throughout the NLR and that they can- 
not, alone, explain all NLR emission because shock signatures 
cannot always be observed. There have been several previous to 
solve the temperature problem from different perspectives (e.g. 



ISinette, Wilson & Storchi-Bergmaniil ( Il996l) ). Here, we compare 
some of the current models with our selected high-Tj object data. 

Line-ratio diagnostic diagrams are commonly used to visual- 
ize and compare models with observations. The ratios of emission 
lines plotted against one other, depending upon the emission lines 
chosen, can show clear relationships in density, metallicity, and the 
ionization mechanism. For example, [O III] /i5007/[O II] /i3726,29 
(hereafter [O IlI]/[0 II]) is theoretically independent of the chemi- 
cal abundance and is a very good indicator of the ionization param- 
eter. [O 1] /i6300/HQ- (hereafter [O I]/Hq') is sensitive to the ion- 
ization parameter. [O III] /15OO7/H0 (hereafter [O llI]/H/3) depends 
on all the parameters but is also frequently used as an indicator 
of the ionization parameter. [N 11] /i6584/HQ' (hereafter [N Ilj/Ha) 
is mostly sensitive to the gaseous abundance. [S 11] /16716,31/Hq' 
(hereafter [S Ilj/Ho-) is sensitive to the density. The ionization pa- 
rameter U is defined as 



U 



4nr-nc 



(3) 



where Qm, is the total number rate of hydrogen ionizing photons 
isotropically emitted by the nucleus, r is the representative distance 
of the emission-line clouds from the central source and n = jih is the 
total hydrogen number density. We calculate the extinction correc- 
tion effects on multiple line ratios: [N llj/Ha, [O llI]/H/3, [S II]/Hq-, 
[O Ij/Ha, [O llI]/[0 II] and the temperature-sensitive [O 111] ratio. 
The [O llI]/[0 II] ratio is the most sensitive to extinction coiTec- 
tion, showing a ~30 per cent decrease after the extinction correc- 
tion, and the temperature-sensitive [O 111] ratio decreases less than 
10 per cent. 

Fig.|7]shows the NLR observations on the line-ratio diagrams. 
Here we select four frequently used line-ratio diagnostic diagrams: 
[N llj/Ha versus [O II1]/H0, [S Ilj/HQ- versus [O llI]AI/3, [O Ij/Hor 
versus [O III]/I^ and [O llI]/[0 11] versus [O llI]/H/3. The obser- 
vational data points are the NLR-dominated objects in A/BSY (i.e. 
sample A/B Seyferts) with > 15 000 K (crosses). The model 
grids are plotted using ite raH which is a IDL tool for emission- 
line ratio analysis (refer to lOroves & Allenllioiol for details). The 
metallicities shown in the side bar are in units of solar metallicity. 

Temperature and ionization are fundamentally deter- 
mined by the ionization source. First, we try both AGN 
(iGroves. Dopita & Sutherland! l2004allbl) and shock ( I Allen et al] 
l2008h model libraries, which are currently available with ITERA. 
The former are photoionization models of the NLRs of AGN, 
while the latter consider fast shock effects. We find that the AGN 
models fit our observations better. Previous observations indicate 
that the NLR clouds are likely to be dusty and clumpy in nature. 
The dusty AGN models consider this effect to produce a physical 
model to explain why and how the AGN NLRs cluster on line 
ratio diagrams: dust dominates the opacity in dusty gas at high U, 
and the gas pressure gradient must match the radiation pressure 
gradient in isobaric systems. Therefore, these assumptions lead 
to a self-regulatory mechanism for the local ionization parameter 
and hence for the emission lines. We have conducted further ex- 
periments and found that the dusty, radiation-pressure dominated, 
photoionized AGN model grids can best reproduce our data when 
assuming a gas density of nH=100 cm"' and a power-law index 
a = - 1 .4. In Fig. [7] we can see that the observed flux ratios are 
generally well described by the dusty AGN model predictions. 



' ITERA, and associated libraries and routines are all available from 
|http://www.brentgroves.net/itera.html| 
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Figure 7. Four panels show selected diagnostic diagrams in which models and observations are compared: (a) the [N II]/Hq' versus [O III]/H/3 diagram; (b) 
the [S II]/Ha versus [O III]/H8 diagram; (c) the [O Il/Hff versus [O IU]/U/3 diagram ; (d) the [O ni]/[0 II] versus [O III]/Hy3 diagram. The observational data 
points (black crosses) are the NLR-dominated objects in A/BSY (i.e. sample A/B Seyfeits) with Tg > 15 000 K. The model grids are plotted using ITERA, 
showing dusty, radiation-pressure dominated photoionized AGN models. We assume electron density of 100 cm"' and the power-law index a = -1.4. The 
metallicity is in unit of solar metallicity. 



However, some remaining problems still need to be addressed, as 
follows. 

(i) Our NLR data present a characteristic range of 10'"' cm"' . 
The fact that the dusty AGN grids fit the data with density = 100 
cm"' better than other density constraints is illustrative. 

(ii) Despite the overall good fitting, the model grids present a 
systematic overestimation by a factor of 2 in the [O III]/[0 II] ver- 
sus [O III]/Hy6 diagram. [O III]/[0 II] is sensitive to extinction cor- 
rection and it decreased by approximately 30 per cent after correc- 
tion. Even when including the uncertainties introduced by extinc- 
tion corrections, the remaining discrepancy is still remarkable. 

(iii) Although emission line ratios support photoionization as 
the dominant ionization mechanism, the high gas temperatures and 
velocities observed within the NLR indicate that shocks might play 
a part in this region (see the discussion in Section 7.2 ). 

(iv) Some of these high-r,, strong [O III]/14363 emission Seyfert 
2 galaxies show low-metallicity (i.e. Z/Z© ~ 1; see the discussion 
in Section 7.3). 

To proceed, we will try to physically model the NLRs and 
probe these problems in a companion paper. 



6 AND SFR/M* 

The electron density and temperature are quantities representing 
the cuiTent statuses of NLRs, while stellar mass (M^,) and star 
formation rate per unit mass (i.e. the specific star formation rate, 
SFR/M,,) provide information about the galaxy formation histo- 
ries. Either stellar mass is conserved or it is a slowly increasing 
qua ntity unaffected by m orphological transformations and mergers 
(e.g. lHopkins et all 20091) . Traditionally, the masses of galaxies are 
estimated by dynamical methods from the kinematics of their stars 
and gas. The specific star formation rate has either explicitly or im- 
plicitly been used in numerous studies of field galaxies at z < 1. 
It defines a useful measure of the rate at which new stars add to 
the assembled mass of a galaxy (i.e. a characteristic time-scale of 
star formation) and it is strongly related to several other important 
physical quantities (e.g. lClark & Oevll2002h . 

We adopt M*, and SFR/ M<. derived using the met hodolo- 
gies oflCauffmann et al. 2003aland lBrinchmann et al.l2004l respec- 
tively. lKauffmann_^^alj (,20^ have constrained the stellar masses 
of galaxies using two stellar absorption-line indices (the 4000-A 
break streng th D„4000 an d the Balmer absorption-line index H<5^). 
[Brinchman n et al.l ( l2004b have obtained the SFR and SFR/M,, by 
two means: for star-forming galaxies, they used the observed Ha 
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Figure 8. Stellar mass versus dust attenuation Ay (top) and specific star formation rate SFR/M* (bottom) in bins of M*. Empty boxes are the median 
values of M-t, Av and SFR/M*, with error bars showing Icr dispersions. From left to right: Seyferts, LINERs, composites and star-forming galaxies. Each type 
of galaxies is divided into sample A (orange), sample B (green) and sample C (blue). Red triangles are the NLR-dominated objects, dark grey open squares 
present the distribution of the sample A objects and grey background indicates the remaining objects (please refer to the online figure for details). 



luminosity and converted this into a SFR, while for composites 
and AGNs, they used a calibration based on D„4000. The sample 
they used was drawn from SDSS DR2 with essentially all aper- 
ture bias removed. N ote that their results are based on a model grid 
jCharlot et al.ll2002h . ensuring a consistent picture for the attenua- 
tion of continuum and line emission photons. 

In this section, we turn our attention to M^, and SFR/M,t and 
we analyse both the relations between the two properties for differ- 
ent classes and their relations with the obtained plasma diagnostic 
results. Then, we compare these physical properties of the NLRs to 
those of the host galaxies. 



6.1 Mj, versus SFR/M,^ 

It is interesting to plot SFR/M* as a function of . Previous stud- 
ies have shown that galaxies generally split into two basic popula- 
tions: concentrated galaxies with low SFR/M, a nd low-mass, less 
concentrated galaxies with hig h SFR/M* (e.g. iKauffmann et al.l 
l2003bl : iBrinchmann et al ] |2004 . However, this is far from the full 
story. ISalim et al. I( l2007l) have constructed a two-dimensional prob- 
ability distribution function in the (log SFR/M», log M*) plane and 
found that different types of emission-line galaxies occupy rela- 
tively distinct portions of the parameter space. Here, we attempt to 
provide more details. 



(i) Fig. [8] (bottom panels) shows that the SFR/My, of the star- 
forming galaxies slowly decreases with increasing M,, from ap- 
proximately 10** to 10" Mq, while those with high M,. show a 
much less pronounced correlation. Interestingly, those 'less active' 
galaxies are indeed active galaxies: Seyferts and LINERs. Most 
composites are located at the junction connecting the AGN and 
star-forming branch. The distri bution pattern agrees well with that 
presented bv lSalim et J]|2007l . 

(ii) As shown in Table [3] objects with high Ay (e.g. sample 
C composites), generally reside in the conjoining region. They 



have a much higher Ay than the Seyferts at constant M*. In 
general, high mass galaxies are older than low-mass galaxies, at 
a fixed stellar mass, studies show that high SFR galaxies yield 
relatively younger populations tha n those with lower SFR (e.g. 
ISchaerer. de Barros &Skliasll2012l) . In other words, galaxies with 
lower SFR/Mj. along with higher M^., on average, have experienced 
acti ve star-forming activ ities and therefore are intrinsically older 
(e.g. lHeavens et alj|2004) . Given that the amount of dust produced 
is proportional to the amount of stellar mass, as shown in the top 
panels in Fig. [8] we speculate that some wide-reaching mechanisms 
(e.g. outflows driven by the AGN) might act to blow out dusty gases 
or consume a significant fraction of the dust content while lowering 
the SFR/M,, for massive galaxies. 

(iii) The selected NLR-dominated objects are plotted in Fig. [8] 
(red solid squares), contrasting with the background host galax- 
ies (grey dots). In the (log M*, log SFR/M^,) parameter planes, 
the NLRs distribute uniformly with the host galaxies, only biasing 
unsurprisingly to higher masses for the composites, because they 
reside in massive active galaxies by nature. To be observed, opti- 
cal SFRs require emission-line luminosities and thus the aperture 
effects on emission-line luminosities can lead to strong biases in 
star-formation rate studies. In the (log M*, Ay) plane, we observe 
that the NLRs-dominated obj ects, on average, p r esent higher Ay 
values, which is also expected. IBrinchmann et al] ( |2004) have pro- 
posed an aperture correction method and have calculated the likeli- 
hood distribution of the SFR for a given set of spectra with global 
and fibre g - r and r - i colours. They have proposed that their cor- 
re ction method is robust, only if > 20 per cent (a criterion argued 
bv lKewlev. lansen & Gelledl2005l) of the total r-band light is sam- 
pled by the fibre; thus, a redshift z > 0.04 is required for SDSS. 
However, as given in Section 5, only systems with ln0 < 1 (i.e. z < 
0.046) are considered to be NLR-dominated. With these considera- 
tions, no precise conclusion can yet be reached, and further specific 
observations might be helpful. 
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Figure 9. Sample distributions on the versus SFR/M^ diagram, where 
red solid squares denote the NLR-dominated objects and grey squares indi- 
cate the remaining objects. The dashed lines indicate log M-,, = 10.48 (i.e. 
M* = 3x10'" Mo yr-l, vertical) and log SFR/M* = -10.30 (i.e. SFR/M* = 
5x10"" yr"', horizontal). The dotted Une presents log SFR/M* = -12. 



(iv) Finally, active galaxies with strong [O III] /14363 emission 
behave unexpectedly. A/BSY and A/BL distribute uniformly with 
CSY/L, and A/BC are located on the high SFR/M^, side with lower 
compared to CC objects (refer to Table |2] for the abbrevia- 
tions). Meanwhile, the situations are quite different in the star- 
forming cases: A/BSF definitely have both much lower M^, and 
higher SFR/M<, compared to CSF objects. As we have previously 
understood, strong [O III] /i4363 emissions in star-forming galax- 
ies imply high Tg and low cooling efficiency, and thus low gaseous 
metallicity; therefore, these objects have higher SFR/M,t than weak 
[O III] /14363 emission star-forming galaxies at constant M^,. How- 
ever, this reasoning will not work in the AGN cases. Taking the 
Seyferts as an example, it seems that the occurrence of strong [O 
III] /14363 emission is just a matter of frequency: there is no signif- 
icant separation between the strong (i.e. S/N > 5) and weaker (i.e. 
1 < S/N < 5) [O III] /14363 Seyferts in the (log M,^, log SFR/MJ 
plane, compared to composites and star-forming galaxies. In ad- 
dition, the strong [O III] /14363 Seyferts distribute in a much wider 
range of SFR/M,^ (i.e. -10.6 < SFR/M,t < -9.2) than the con-espond- 
ing star-forming ones (i.e. -9 < SFR/M^. < -8.2). Assuming that 
the classifications of nan'ow emission-line galaxies are correct, we 
believe that such significantly different distribution patterns reveal 
some intrinsic physics related to [O III] /14363 generating mecha- 
nisms and the properties of the central energy sources in galaxies. 
Some studies have shown that, in general, the ionization potential 
of the e mission region is positively correlated w i th the [O III] A436 3 
flux (e.g. lNagao. Muravama & Taniguchil200ll : l\^ona et alj2012l) . 
Therefore, we refer to the results described above as the 'evolution- 
ary pattern of AGN with high ionization potential' . The key point is 
regarding the nature of the high-energy processes that excite the [O 
III] corona l lines in active galaxies. For example, the results from 
lAird et al] ( l2012h demonstrate that the same physical processes reg- 
ulate AGN activity in all galaxies in the range (i.e. 9.5 < log 
M^, < 12) and could most likely be provoked by the energetically 
central nucleus. 



Table 9. Summary of SFR/M* in seven bins. 
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6.2 Galaxy formation timescale sequence 



It is now commonly accepted that the star formation history of a 
given galaxy depends strongly on its mass. Higher values of the 
specific SFR suggest that a larger fraction of stars was formed 
recently. The specific SFR is often called the galaxy present-day 
build-up time-scale, while is the product of the overall galactic 
star-forming processes. We compare the mean SFR/M^, of the four 
classes in seven bins. Table|9]clearly demonstrates that SFR/M/. 
> SFR/Mi K < SFR/Mc < SFR/Ms/r at constant M,. Therefore, we 
consider these two clues combined as an indicator of galaxy forma- 
tion history: it is generally true that Yt > Ysv > Yc > Ysf, where Y is 
present-day star-formation time-scale. However, this should not be 
confused with an actual galaxy age, because the specific SFR tells 
us only how long it would have taken to build a galaxy, assuming it 
had a current SFR throughout its lifetime. 

In Fig. [21 we plot four classes of emission-line galaxies: 
Seyferts, LINERs, composites and star-forming in the (log M^, log 
SFR/M*) plane. We define an empirical base point on the logscale 
M,, versus SFR/M* diagram, [-10.30, 10.48], which corresponds 
to = 3 X 10'" Mo and SFR/M* = 5 x 10"" yr"' (marked in 
Fig. |9](. The latter threshold is a 'by-eye' division, because most 
star-forming galaxies perfectly distribute above the vertical dashed 
line. lKauffmann et al. ( 2004,) have found that the transition from a 
young to old stellar population takes place at a characteristic M^. of 
3 x 10'" Mo, the logscale value of which is ~ 10.48. For galaxies 
with , above this threshold, the distribution of sizes and concen- 
trations at fixed stellar mass is independent of their local environ- 
ments, and these galaxies have high stellar masses, high concentra- 
tions, high surface densities, little on-going star formation and red 
colours. 

Now we use the base point described above as the vertex to 
determine the angle, which is labelled ipsFR/Mt' that each galaxy 
in our sample makes with the x-axis in the (log M*, log SFR/M^.) 
plane. Fig.llOlshows the normalized histograms of this angle corre- 
sponding to different galaxy types. Notice that, for accuracy, those 
with log-scaled specific SFR exceeding -12 are not included (i.e. 
objects that lie below the dotted lines in Fig[9]are not counted be- 
cause SFR/M^ estimate s trend to be inaccurate past ~10"'- yr"'; 
iBrinchmann et al.l I2OO4I) . We can see that the histograms of the 
Seyferts and LINERs look quite similar and both peak at <psFR/Mt 
~ 300°, while composites and star-forming galaxies specify very 
different distributions: the former peaks at approximately 90° and 
the latter at approximately 140°. 



6.3 Current status versus formation history 

Utilizing previously obtained and T^, which indicate the current 
statuses of galaxies, along with the dust attenuation Ay, finally we 
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Figure 10. Nomialized histograms of (psFR/M-,,, which is the angle relative to the x-axis in the versus SFR/M^ diagram, using the empirical point [-10.30, 
10.48] as the vertex (see Fig. |9). From top to bottom, we show Seyferts (A/B/CSY; refer to Table |2]for the abbreviations), LlNERs (A/B/CL), composites 
(A/B/CC) and star-forming galaxies (A/B/CSF). Blue dotted histograms show objects with 1 < S/N < 3 for [0 III] /14363, green dot-dashed ones indicate those 
with 3 < S/N < 5 and orange solid (shaded) histograms are those with S/N > 5. Note that objects with log SFR < -12 are not included. 



analyse their relations with galaxy formation history, indicated by 
and SFR/M^.. Fig.llllillustrates our sample distributions in five 
property pairs. We notice the following. 

Dust attenuation. The amount of dust content shows a pro- 
nounced correlation with the amount of stellar mass (i.e. massive 
galaxies contain more dust). As a function of SFR/M,t, Ay shows 
a positive correlation at low SFR/M^, and then becomes negative. 
Objects with log SFR/M*, ~ -10 have the highest A v, and the NLR- 
dominated objects (red triangles) have a wide range of Ay and show 
both high Mj, and low SFR/M* , compared to the host galaxies (grey 
dots). 

Electron density. As a function of , n^ also shows a positive 
correlation and reaches the highest value at the massive end. As 
a function of SFR/M,., a significant threshold is clearly revealed. 
A large fraction of objects with log SFR/Mj, < -10 have n^ > 100 
cm"\ while those with log SFR/M* > -10 show a much wider n^ 
distribution. Most of the NLRs naturally lie in high regions. 

Electron temperature. Compared to both Ay and n^, displays 
much less pronounced coiTelations with both M,, and SFR/M^. The 
existence of an unexpected high-r,, branch is quite remarkable. The 
branch appears at approximately 9.5 < log M,. < 1 1 and -10.4 < log 
SFR/M,, < -9.4. The NLRs as a group show higher than the host 
galaxies. 

A significant T,, 'pumping-up' is shown in Fig.[TT] As listed in 
Table. [To] while having similar redshifts, objects with T,, > 2x10* 
K show much lower [O III] /15007. Considering that increases 
rapidly with decreasing R[0 III] (i.e. /[O III] /i5007/.i4363), we 
infer that the plume towards high is most likely the sign of an 
inaccurate estimate at the low R[0 111] end (Fig.ll2b. Moreover, 
we can see in Fig.|TT]that objects from sample A and sample B form 
a bimodal distribution in M-,, and SFR/M,,, which is not seen in 
sample C. Given that [O III] /i4363 correlates with galaxy activity, 
we argue that there exists a distinct low-activity population in active 
galaxies with 9.5 < log M,^ < 11 and -10.4 < log SFR/M* < -9.4, 
which might indicate some evolutionary stage. 



Table 10. Summaiy of high- objects (> 2x10* K) and low-r^ ob- 
jects (the rest) from samples A and B (median values). Note that [O 
11I]63 and [O llI]o7 denote the pre-extinction correction fluxes of [O 
111] /14363 and 5007, respectively (in units of 10"'^ erg/s/cm"). 
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7 DISCUSSION 

7.1 Shock effects on n^ and 

Stellar photoionization is not the only process that leads to the for- 
mation of emis sion lines. In fact, ionization can al s o be p roduced 
in shocks (e.g. lOopita, Groves & Sutherlandlll995L Il996l) . To ex- 
plain the observed emission line ratios, the possible contribution 
of shock-heating to photoionization must be examined. The cool- 
ing radiation from shocks produces high-energy photons that con- 
tribute to ionization. Kinetic energy is transformed into thermal en- 
ergy, and the energetic thermal electrons can ionize the gas. Shocks 
are ubiquitous in galaxies. For example, shocks in AGNs can be 
caused by jets, winds or supersonic turbulence. Here, we find sev- 
eral emission line spectral indicators that can help trace the effects 
of shocks. 

(i) The density calculations show that the strong [O 111] /i4363 
AGN density is 100-150 cm""* higher than those that show no [O 
III] (see Section 4). Shocks generate compression, and therefore 
high gas densities. While temperature-sensitive line ratios change 
with the metallicity of the gas, they are also strongly affected by 
the ionization state of the gas. Therefore, shocks and the varying 
contribution of star formation and AGNs to the ionization become 
important. 

(ii) Table|7](see column 12) lists objects with Tg > 20 000 K in 
each subsample, and we have shown in Section 5.3 that the dusty 
AGN model grids can fit some of these objects well. Shock heating 
produces very high temperatures (of the order of millions of K), 
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Figure 11. Distribution of galaxy properties. The diagonal panels show the nonnalized distribution of five properties independently: M*, SFR/M*, Ay, Tg and 
«(, (grey solid ones, sample A and sample B; lightgrey dot-dashed histograms, sample C. Objects from samples A and B present a bimodal distribution in M,^, 
SFR/M-), is apparent, which is not seen in sample C. The off-diagonal panels show the bivariate distribution of each pair of properties, revealing the complex 
relationships among them. Red triangles are the NLR-dominated objects, blue open boxes mark out LINERs, while grey dots present objects from samples A 
andB. 



leading to collisional ionization, and the production of highly ion- 
ized species. In fact, as we can see from Table|7]and Fig.[TT] high 
Tg species are more easily found in LINERs and composites with 
strong [O III] /14363 emission with certain SFR/M* values. Thus, 
this indicates high-energy processes in their emission-line regions 
and the significant contribution of shocks evoked by the active nu- 
cleus, and even hot post-asy mptotic giant branch stars and white 
dwarfs must be important (e.g.|Stasinska et al. 2008). 

iGroves. Dopita & Sutherland! ( l2006ah have reported that the 
observed NLR spectra reveal much flatter mid-IR spectral slopes 
than those obtained in most NLR models. This suggests that shocks 
are a possible explanation for the difference between the model and 



observed slope. We confirm that it is both sensible and necessary to 
consider the combination of photoionization and shock. Additional 
multi-wavelength (especially infrared) observations could be very 
helpful in achieving this goal. 

7.2 Strong [O 111] narrow-line AGN: low-metallicity 
candidates? 

"The physics of the NLR is simpler than that of the broad-line re- 
gion; the physical conditions are in fact similar to H II regions, so 
that many of the metallicity-sensitive emission line ratios that are 
routinely applied to star-forming galaxies are also good diagnos- 
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Figure 12. Top: [O III] A5007 versus [O III] i4363 (pre-extinction correc- 
tion fluxes in units of 10"'^ erg/s/cm-). Red squai'es are objects with > 
2x10^ K and grey crosses present those with lower T/,. Bottom: normalized 
histograms of [O III] /15007 fluxes for the above two groups. Note that only 
objects from samples A and B are plotted. 



tics in thie NLRs. IStorchi-Bergmann et al.l (Il998h have calibrated 
several of tlrese emission line ratios for nearby AGN, using H II 
region determined metallicities and photoionisation models. How- 
ever, the similarities between the NLRs and the H 11 regions in a 
galaxy lead to ambiguities when estimating abundances. Unless 
spatially resolved observations can be obtained, it is often difficult 
to disentangle the contributions to the emission-line spectrum from 
the NLR and star-fo rming regions within the galaxy", cited from 
lGrovesetal] ( l2006d) . 

Low-metallicity AGNs appear to be very rare objects, at least 
locally. Currently, emission-line-based estimates of AGN metal- 
licities at both high and low redshifts indicate that AGNs have 
predom inantly solar-to-supersolar metallicities (e.g. Ig roves et al.l 
l2006bl) . This phenomenon is expected only if AGNs are old galax- 
ies, are experiencing more supernova explosion events or even 
mergers, and are sitting on the later part of galaxy evolutionary 
sequence. From a theoretical point of view, when the total metal- 
licity of a photoionized nebula decreases, there are several ef- 
fects that cause changes in the final emission-line spec trum (e.g. 
iDopita & Sutherl and 2003; lOsterbrock & Ferlandll2006l) . For ex- 
ample, as metallicity decreases, the temperature of the nebula con- 
versely increases. This is a result of the decrease in the efficient 
cooling metal emission lines and the availability of more high - 
energy photons to ionize the hydrogen ("Sutherland & DoDital2003h . 

Using the SDSS DR4 data, Groves et al. (2006b) have found 
40 low-metallicity AGN candidates selected out of the 23 000 
Seyfert 2 galaxies. They expected that these candidates with masses 
below 10'" Mq could have metallicities approximately half those 
of typical AGN. Coincidentally, many of the candidates they found 
show the [O III] /14363 line in their spectra. They have suggested 
that the use of a low-mass selection criterion supported by the ob- 
servations reveals a correlation between mass- and metal-sensitive 
line ratios, and that the existence of [O 111] /14363 indicates higher 
gas temperatures in these objects and thus low metallicities. If 
their results are correct, we suggest that our strong [O 111] /14363 
emission AGNs selected from SDSS DR7 similarly contain low- 



metallicity AGN candidates (as shown Fig. O. Using our refined 
ranges of density and temperature as input parameters for our NLR 
modelling, which will be the topic of a companion paper, we aim to 
reproduce the observational evidence and, as a definite challenge, 
obtain metallicity calibrations for active galaxies in order to further 
check these candidates. 



8 SUMMARY 

We have presented a statistical study based on spectroscopy of a 
large-sample from SDSS DR7. Our sample comprises 15 019 ob- 
jects with a median redshift ~ 0.08 and is divided based on [O 
II1]/14363 detection quality: 1100 in sample A (i.e. S/N > 5), 1411 
in sample B (i.e. 3 < S/N < 5), and 13 212 in sample C (i.e. 1 < 
S/N < 3). First, we determined electron density and electron tem- 
perature through /[S II] /16716 //16731 and /[O III] /15007//14363, 
respectively. To focus on the NLRs, we required that FWHMBaimer 
> 300 km s"' and the physical aperture size Incp < 1 kpc to select 
the NLR-dominated objects within our sample. Then, we analysed 
a typical range of rig and of these objects and compared them 
with those of the host galaxies. Furthermore, we have utilized the 
plasma diagnostic results along with and SFR/M<- to study the 
relations between these physical properties for four different classes 
of emission-line galaxies. Our main results are listed as follows. 

(i) The typical range of density in NLRs of AGNs is 10^"^ cm"'. 

(ii) The typical range of temperatures in NLRs of AGNs is 1.0- 
2.0x10* K for Seyferts; for LINERs and composites, this range 
could be higher and wider. 

(iii) Transitions of both iig and from the NLRs to the discs are 
revealed (i.e. both and decrease outwards). 

(iv) Two sequences are proposed: Tuner ^ T,„„,p„,i,g > Tgey/ert 

^ Tslar- fanning and llLlNER <^ f^Seyfer! ^ f^composire ^star-forming • 

(v) The median value of «(, is ~ 150 cm"' higher in both Seyferts 
and LINERs with strong [O III] /14363 emission than in those with 
weak [O III] emission, most likely indicating some effect of shocks. 

(vi) LINERs and composites show the median to be approx- 
imately 2.0x10'' K, which is twice as much as that of star-forming 
galaxies and far too high to be explained by only stellar photoion- 
ization. 

(vii) Some strong [O III] /i4363 emission Seyfert 2 galaxies with 
> 15 000 K can be fitted with dusty AGN model grids at low- 
metallicity (i.e. Z/Zq ~ 1). 

(viii) LINERs with strong [O III] /14363 emission simultane- 
ously contain more dust and indicate the coexistence of old stellar 
populations and some specific high-energy process. 

(ix) Objects with high Ay reside in a certain region of the 
M^, versus SFR/M^, diagram, suggesting that some wide-reaching 
mechanisms blow out or consume a significant fraction of dust con- 
tents, and meanwhile lower SFR/M* for the massive galaxies. 

(x) Considering characteristic present-day star-formation time- 
scale Y, we suggest the following relationship: Yi > Ksi' > Kc > 

(xi) Seyferts with strong [O III] /14363 emissions lie in the high 
M^, region and distribute uniformly with those showing weaker [O 
III] /i4363 emissions (the so-called evolutionary pattern of AGN 
with high ionization potential). The central nucleus might be re- 
sponsible for the high ionization potential in Seyferts. 

(xii) The NLR-dominated objects have higher M^, and lower 
SFR/M^. than the host galaxies, most likely suggesting that the 
NLRs reside in old and massive active galaxies. 
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In conclusion, we find several pieces of evidence pointing to 
additional source(s) of ionization other than the fact that photoion- 
ization by hot stars must be working to produce the highly ionized 
species such as the strong [O III] /14363 line in NLRs of AGNs. 
Previous observations have shown correlations between jets and 
NLRs. The high gas temperatures and velocity dispersions present 
in our sample imply that the combination of shock and radiation 
heating might work as well. The best-known choice of ionization 
source is the energetic central nucleus, and (weak) shocks might be 
generated by feedback from AGNs in mergers or by the inflation of 
cavities by the central AGN. With our refined ranges of density and 
temperature as input parameters for the AGN NLR modelling, in a 
companion paper, we intend to reproduce the observational spectra 
and, furthermore, to obtain metallicity calibrations for active galax- 
ies. 
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